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SUMMJty 


A  nMthod  la  davelopad  to  find  tha  atraaaaa  and  atralna  In  an  Incoapraa- 
alble  vlaco-alaatlc  hollow  cyllndar  with  annihilating  innar  radiua  contained 
by  an  alaatlc  caaa  and  aubjact  to  Internal  praaaura  under  the  aasumptlon  of 
a  atate  of  plane  atraln. 

Streaaea  and  atralna  are  conputad  foe  a  natarlal  with  davlatorlc  atreaa- 
atraln  relatione  characterlatlc  of  a  atandard  aolld.  The  nunarlcal  computa¬ 
tion  la  carried  out  with  the  aid  of  an  I.B.M.  digital  coag»uter  1620  and  la 
intended  to  llluatrate  tha  effacta  of  tha  thieknaaa  of  tha  cyllndar,  of  tha 
rata  of  Incraaaa  of  tha  Intamal  praaaura  and  of  tha  atrangth  of  tha  raln- 
forceoMnt  provldad  by  the  elaatlc  ahell. 
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1.  Introduction 


The  problem  of  evaluation  of  stresses  and  strains  In  a  solid  propellant 
with  a  burning  Inner  surface  under  Internal  pressure  has  not  yet  been  solved 
In  a  general  way  because  of  the  evident  difficulty  of  the  treatment  of  the 
boundary  condition  at  the  moving  Inner  surface. 

It  is  the  purpose  of  the  pre;>ant  Investigation  to  find  a  solution  of 
this  problem  by  evaluating  stresses  and  strains  In  an  elastically  case-bonded 
linear  vlsco-alastlc  hollow  cylinder  with  annihilating  inner  s;irfaca  subject 
to  Internal  pressure. 

For  simplicity,  incompressibility  of  the  material  and  a  state  of  plain 
strain  are  aasumsd;  hence 

e,  *  •  u/r 

where  V  -r  and  e^  are  raspaetively  the  radial  and  tangantlal  compaasnts 
of  strain  and  devlatoric  strain,  and  u  ia  tha  radial  displacosMnt. 

Eq.  (1)  is  satisfied  by  u  ■  k(e}/r  where  k(t)  is  s  function  of  time  t 
only.  Uanca, 

C^--kCi)/r*  g*-  kCtVr*  (2) 

When  s^j  and  a^j  raapactivaly  denote  the  components  of  devlatoric  streas 
and  strain,  the  stress-strain  relations  for  incompressible  linear  visco-elaa- 
tic  materials  are  of  the  well-haown  form 
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wiMr*  P  and  Cl  arc  llnaar  diffarantlal  oparatora  with  raapact  to  tlaa  t  or, 
in  taraa  of  Laplaea  traaaforaa  with  Initially  aaro  eondition, 

(5') 

whara  p  danotaa  tha  tranaform  paraMtar. 

Tha  ofuation  of  aquilibriua 

•^  + o  or  4-  ~  ^  ■  O  (^) 

whan  P  oparataa  on  both  sidaa  with  aaaiwiad  intacehangaability  batwaan  dif- 
farantiationa  with  raapact  to  tiaa  t  and  radiua  r  baeonaa 

•  |ClCi^-  CKap]/r 

after  tha  atraaa-atrala  relation  Eq.  (3)  baa  bean  introduced. 

Subatitutlag  Bq.  (2)  into  the  abova  aqpaation  givaa 

-J:{  ipCO]  • 

fron  which  IP(o^)  ia  obtained  by  intagration  with  raapact  to  r: 

~T  *  cc«  (5) 

where  C(t)  ia  a  function  of  t  only. 

Tha  aquilibriua  of  tha  alaatic  aball  raquixae  tha  relation 

•-rlr-b  • '■b'«c 

batwaan  tha  ahell  atroaa  0^  and  tha  radial  atraaa  of  tha  eyliadac 
at  tha  interface  batwaan  tha  cylinder  and  tha  ahall,  idiila  o^  ia  in  turn 


I 


related  to  the  strain  of  the  cylinder  this  interface  by  the  rela¬ 

tion 

*  I  ^ >^0*  ®’e^r*b 

since  s^]  .  is  identical  to  the  shell  strain  e  :  h  is  the  thicknees  of  the 

e  rab  c 

shelly  b  the  outer  radius  of  the  hollow  cylinder  and  E  and  k*  are  Young's 

c  o 

modulus  and  Poisson  ratio  of  the  shell  raspaetivaly. 

Therefore  the  first  boundary  condition  at  r  ■  b  has  the  form 

®r]r.W  “  “ 
or 

1PMr.b*  («•) 


where  the  second  of  Bqs.  (2)  has  been  introduced  for  end 


e: 


Efc 

I  -  K 


IT  ♦ 


(7) 


C(t)  is  obtained  from  Eqs.  (^)  and  (6'): 

ca)  -trUlkC^  -  iPjWCft} 


Hence 


WW  (tf-  .  (8) 


However,  as  pointed  out  by  B.  H.  Ine,  R.  M*  Radok  and  H.  B.  Woodward  [1], 
it  is  not  possible  to  operate  with  IP  on  the  second  boundary  condition  at 
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the  Inner  surface 


-•’rlr.eaVN*  (9) 

and  substitute  It  Into  Eq.  (?)  as  in  the  first  boundary  coBdittan,  since 
£q.  (J)  is  essentially  a  relation  batwaan  strass  and  strain  aasoelatad  with 
a  fixed  material  particle,  while  at  each  Instant  ^  >4*  (9)  rapra- 

sents  the  radial  straas  of  a  dlffarant  aatarlal  particle;  a(t)  la  a  mono- 
tonlcally  Increasing  function  of  tlrna  rapraaantlng  the  lanar  radius  of  tha 
hollow  cylinder  until 

(10) 


is  reached  when  tha  annlhllatlao  la  coiBlatad  at  t  ■  t  . 

o 

2.  General  aitpy^aian^  for  gtraaaas  and  Strains 

The  difficulty  Is  raaovad  however  If  the  straas-straln  relation  la  ax* 
pressed  in  tha  fora 

(U) 

instead  of  Bq.  (}),  where  ^  la  a  linear  operator  In  derivatives  and  Integrals 
with  respect  to  tlas  only,  or  In  terns  of  Its  Lsplare  transform  under  aero 
initial  conditions: 

(11*) 

The  substitution  of  Iq.  (2)  Into  Iq.  (11)  gives 
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The  radial  streaa  is  then  obtained  from  Eq.  {k): 

<r^  «  +  Ort') 

where  D(t)  is  a  function  of  t  only. 

Uae  of  the  firat  boundary  condition  Eq.  (6)  determinea  0(t): 

-  e: 

Hence 

•  (-t* ' -  K ^  (>2) 

The  unknown  function  k(t)  ia  now  obtained  from  the  lntegro>dif£erantiel 
equation 

(13) 

tfhich  la  the  reault  of  the  iqiplieation  of  the  aeeond  boundary  condition 
Eq.  (9)  to  Iq.  (12). 

B.  H.  Lae,  R.  H.  Radok  and  H.  B.  Woodward  [1]  aolved  a  apecial  case  of 
the  problan  in  which  the  aecbanical  aodel  of  the  laaterial  was  aaauned  to  be 
a  Kelvin  body.  This  implies  a  stress -a train  relation  of  the  form  Eq.  (11) 
with  _ 

^  »  2G  TF 

or  of  the  fora  of  Eq.  (11*)  with 

2G|♦^^^p; 

G  denotes  the  shear  aodulus  and  Ij  the  coefficient  of  viscoaity  in  shear  as 
shown  in  Fig.  1. 

In  aore  generality  linear  visco-elastic  aster ial  can  be  represented  [2] 
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by  n  Maxwell  elemanta  coupled  In  parallel  with  diacrata  ralaxaelon  tlaaa 
\  (Fig.  2)  and  noraaliaad  diacrata  ralaxation  apactrvm  F^(t^)  ■  G^(t^)/G 

(Fig.  })  or  by  the  limit  of  the  same  modal  whan  n  approachaa  infinity  with 
continuous  ralaxation  time  t  and  normaliaad  continuous  ralaxation  spectrum 
£('^}  ■  8(y}/G  where  the  unralaxad  shear  aK>dulus  G  ■  G  ■  g(T)dT. 

It  can  than  be  shown  [3I  that  the  strass  deviation  is  related  to  the  strain 
deviation  in  the  form 

«  lGil«q  ♦/  (Ih) 

at  in  tha  form  of  Bq.  (11')  with 

igCp)  -  . 

The  ralaxation-rata-function  l^(t)  ■  dl^(t)/dt  (this  notation  for  tha  time 
derivative  is  used  henceforth)  in  Bq.  (Ih)  is  the  time  derivative  of  the  re¬ 
laxation-function  ^(t)  idiich  is  defined  so  as  to  produce  the  strees  reaponse 
to  the  unit  step  strain  input  •£j(t)  <■  *0  ■  const,  in  the  form  •j^j(t)  ■  SOaJ^t). 
^(t)  can  also  be  obtained  from  a  knowledge  of  tha  relaxation  epeetron: 

(15) 

*.•1 

for  the  discrete  relaxation  spectrum  and 

1  (16) 

j 

for  the  continuous  relaxation  spectrum.  Bq.  (12)  can  now  be  written  in  the 
form 

r,.- (12*) 

o  * 

I 

j 

i 

I  6 


I 


The  determining  eqtiation  Eq.  (I3)  for  k(t)  than  bacoMS 

which  can  be  reduced  to  the  form  of  Voltarra's  integral  aquation: 


(13*) 


k(t^ 


*  1  .  y  ’ 


i»a-€Dk(ida6' 


where  u  >  B"/aG. 
c 

The  aolutlon  k(t)  of  Eq.  (I7)  fumlahas  the  stralna  with  the  aid  of 
Eq.  (2).  The  radial  atraas  o^  la  obtained  from  Eq.  (12*),  which  iavolvaa 
the  Integration  of  k(t).  However,  the  Identity  from  gq.  (13* ) 

2erCW4^+^Va-fl^k(»a*]-[^  (18) 


can  be  uaad  to  avoid  the  evaluation  of  thla  latagral  writing  the  radial  atreaa 


o^  la  the  fora 


I  I 


«r  •  t-§-w»  -  pce]  -  -^wte  ft .  o 

"C'’”  a*W 


(19) 


Tha  tangential  atraea  9^  la  obtained  from  the  aquation  of  equilibrium 


Eq.  (4): 


(20) 
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when  t  approachea  the  right  hand  aide  of  Eq.  (18)  should  be  Interpreted 
os  the  limit: 

B* 

*  lGr[k(jt;V/  U»"  %  **  r~ —  (21) 

• 

since  “  **  ti?to®c  “  tJStoP^^)  ***■  (^3'). 

Hence,  as  t  approachea  the  atreaaea  are  given  by 

(22) 

V"^  (23) 


since  r  «  b. 

V-^Wq 

3*  Streaaea  and  Strnina  for  Mediua  with  Standard  Solid  Streaa-Straia  Rela¬ 
tion 

When  the  a trass -a train  relation  in  shear  is  rapresantad  by  a  standard 
solid  (Fig.  4),  the  relaxation  tiae  is  considered  infinite  tihila  Tg  ■ 

The  relsxstion-rate-futtction  P(t)  is  obtained  from  Kq.  (1^) 


Hlh 


(15') 


where  sca  >  ®1^(®1  + 

The  substitution  of  Eq.  (1^')  into  Bq.  (17)  produces 


w« 


Jm' 


(17') 
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i)h«r«  h(t)  ■ 

Both  sldoa  of  Bq.  (17* )  or*  now  dlvldod  by  tha  eoofflcimt  of  tho  Into- 
gral  on  tha  laft  hand  alda  and  diffarantlatad  onea  with  raapact  to  tiiaa  lo 
aa  to  produca  tha  diffarantlal  aquation 


whara 


(zk) 


(25) 


and 


(26) 


Hlth  aa 


Intagratiag  factor  axp  [-^  fT^^'  ***  ^  latagratad; 

W«  . 


and  tharafora 


(27) 


It  la  aaauaad  that  tha  aonlhllatlnf  rata  of  tha  Inaar  aurfaca  la  govacnad 
by  tha  xalatlou 

i-Kt  (aB) 
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with 


where  is  the  Initial  Inner  radius  (at  t  ■  O)  and  m  b/a^.  a(t)  Is 
plotted  against  t/t^  In  Fig.  ^  with  Ji^  (■  1.^,  2.0  and  3)  as  paraaster. 

with  the  form  of  a(t)  given  In  Eq.  (26)^  the  Integral  ^  £^)  ^ 

evaluated  so  that 


t±=. 


0 


where 


If  the  Internal  preesure  p(t)  is  assuaed  to  be  of  the  fora 

•  p.O  - 


Bq.  (27)  becoaee 
/.  -t/r, 


ZG(A-  ^  JUe  “T*A'  X;  «  / 

which  can  be  transfocaed  by  Introducing  the  non-diawislonel  quantities 
0  -  P,,/(aG),  de'  »  Tg/T^  end  ;i .  t^/Xg  Into 

^  f.'-' 


-ei'AT^ 


I 


«her« 


Th«  Inecgral  ^^(t)  can  b«  nuMrically  ifvaliuiettd.  Th*  atralna  arc 


*r* 


kft>  Wft^ 


(33) 


whara  R  ■  r/a  .  whlla  eha  atrasaaa 
o 


'Vl  -  «'•*'') f(0  i 


>  (39) 


Eq.  (39)  haa  baaa  afctataad  eonaidarlat  that  by  I'Hoapltal'a  rula 

4.  Nu|SEi££L.SE9BiS 

Tha  pravloualy  dafinad  oaB»ill«anaianal  qaantltlaa  ahieb  eharaatarlaa 
tlM  problaa  ara  vaviawad: 

/>  ■  b/a^  rapraaaata  tha  acl«iaal  thtcimaaa  af  tha  cyltadar. 

O  9 
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X  B  t^/Tg  shows  how  nuch  time  is  required  for  Che  toCel  snnihiletion 

in  terms  of  the  relexstion  tiise  of  Che  meteriel, 

oc'  ■  '^2^^o  ^  measure  of  the  build-up  of  internal  pressure  caaq>ared 

to  the  relaxation  time  of  the  material. 

Amc'  B  c  /t  therefore  indicates  how  rapidly  the  internal  pressure  is 
o  o 

built  up  to  Che  msximum  in  comparison  with  the  total  annihi¬ 
lation  time  t  . 

o 

eta  -f  Gg)  is  releted  to  Che  stress  relaxation  of  Che  standard 

solid  material  in  such  a  way  chat  the  penaenent  screes  after 

completion  of  the  stress  relaxation  under  constant  strain  e^ 

is  s  a  aCe  at  idiile  the  initial  stress  is  s.  a  2Ge  . 

00  o  o  o 

a  a  p^/2G  relates  the  magnitude  of  the  mexlisum  internal  pressure  Co  the 

unrelaxed  shear  modulus  of  the  material. 

E**  20  h 

^  a  ^  a  ^1^^  -  ^  )  *  b  Strength  of  the  reiafoccament  due 

Co  Che  elastic  shell  where  G.  sad  are  the  shear  modulus 

e  c 

and  Poisson  ratio  of  the  elastic  shell  respectively. 

The  assigiaMot  of  nuswrlcsl  values  for  these  quantities  is  listed  in 
Table  1  producing  six  different  cases  for  which  the  numerical  cosiputations 
have  been  carried  out  using  a  1620  I.B.M.  digital  computer. 

Cases  I-lIl  are  for  coavarison  of  the  effect  of  the  origiasl  thickness 
of  Che  cylinder.  The  coovsrison  of  Cases  1,  IV  sad  V  shows  Che  effect  of 
either  the  rate  of  the  application  of  the  internal  pressure  when  t^  is  iden¬ 
tical  in  these  three  cases  (since  than  in  Eq.  (30)  is  equal  to  t^no, 

and  t^  respectively)  or  the  tism  t^  of  total  annihilation  of  Che  cylin¬ 
der  when  -  Tg  is  identical  (since  than  c^  .  lOr^  yt^  and  respectively). 


Finally  Cm«s  I  and  VI  llluatrata  tha  rainf arcing  action  of  cha  olaatlc  ahall. 
In  tha  latter  tha  alaatlc  ahall  la  ao  atrong  that  It  glvoa  rlaa  to  nagatlva 
tangential  atraaaaa  In  tha  cyllndar.  Flga.  6(a),  (b)  and  (c)  raapactlvaly 
ahow  tha  affact  of  the  original  thlcknaaa  of  the  cyllndar,  tha  rata  of  the 
load  application  and  of  tha  rigidity  of  tha  elastic  shall  on  tha  radial  strass 
o^  while  Figs.  7  and  8  show  those  affects  on  tha  tangential  strass  and  strain 
respectively.  Fig.  9  ahoNB  tha  shall  stresses  vlth  h/b  ■  100  for  tha 

six  cases  as  functions  of  tisM  t/t^.  Finally  tha  space  distributions  of  tha 
radial  stress  for  these  six  cases  ara  also  plotted  at  various  tliMS  t/t^  In 
Fig.  10,  while  In  Figs.  11  and  12  tha  spaoe  distributions  of  the  tangential 
stress  and  strain  are  shown.  In  Figs.  10-12,  vertical  straight  lines  corre¬ 
spond  to  the  positions  of  the  ineer  surface  at  the  specified  tiaa. 

TABU  1 


Designation 

/o 

X 

I 

1.3 

10 

.5 

«  ■  *5 

11 

2 

10 

.5 

ec'.  1.0 

III 

3 

10 

.5 

0  ■  .1 

IV 

1.5 

5 

.5 

for  all 

V 

1.5 

1 

.5 

cases 

VI 

1.5 

10 

If 

13 
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Figure  7*  The  Tangential  Streaa  at  the  Inner  Surface 
and  the  Interface  aa  a  Function  of  Time  t/t. 
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20 
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